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Abstract 

An approximate inverse scattering method ^ has been used to construct separable 
potentials with the Laguerre form factors. As an application, we invert the phase shifts 
of proton-proton in the ^5*0 and F2 channels and neutron-proton in the ^Si Di 

channel elastic scattering. In the latter case the deuteron wave function of a realistic np 
potential was used as input. 



1 Introduction 

Low-rank separable potentials greatly simplify many-body computations. The main sources 
of separable nucleon-nucleon {NN) potentials are the inverse scattering methods (see |1| and 
references therein). On the other hand the model effective two-particles interactions (see e.g. 
0) constructed with use of a realistic A^A^ potential (e.g. |^ ) as input are also of separable 
form. 

Recently a novel method was proposed for the treating of three-body problems including 
Coulomb forces [Q. In the framework of this approach the potential operator of the short- 
range part of full interaction is expanded with use of Coulomb-Sturmian (Laguerre) functions. 
Whereas Coulomb interaction is kept in the Green's operator involved in Lippman-Schwinger 
equation. Instead of separable potentials (of relatively high rank) resulting from a realistic 
A^A^ interaction, the potentials obtained by means of inverse scattering technique ^ within 
so-called J-matrix method can be applied in the three-body calculations [Q. 

The J-matrix method 0, is equivalent to the potential separable expansion approach 
|T0| , p!T| ; here too the potential V is approximated by its projection onto the finite subspace 
spanned by the first A^ (Laguerre or oscillator ) basis functions. However in the J-matrix method 
the Schrodinger equation rather than Lippman-Schwinger one is solved in the square integrable 
basis representation. In the J-matrix approach an auxiliary function have been introduced 



resulting from discrete representation of the full Green's operator |T2[. This function (or matrix 
in the multichannel scattering case), so-called P-matrix, plays in our inversion scheme a part 
similar to that of the i?-matrix in the potential reconstruction in the framework of finite- 



difference analogue of i?-matrix scattering theory |jT3[. In particular, information requisite for 



the potential construction can also be obtained from the poles and residues of the P-matrix. 

As an application of our method, we invert phase shifts of the proton-proton and neutron- 
proton elastic scattering. In the A^A^ potential construction an attempt have been made to 



account for the deuteron wave function using the phase-equivalent transformation |T^, ^ . 

In Sect. 2 we outline the J-matrix formalism for the solution of multichannel scattering 
problem. In Sect. 3 we briefly present the inversion procedure in the Laguerre form factors 
case, and describe the taking account for deuteron wave function technique within our method. 
We then apply the method to the A^A^ elastic scattering in Sect. 4. In Sect. 5 we summarize 
our conclusions. 



2 Elements of the J-matrix formalism 



The formalism is spelled out in detail in for treating the problem of multichannel scattering 
from the potentials of the form 



v= 

o, /3 = 1,2 

where the partial wave potentials ^ are given by the expansion 



(1) 



2 Na-l^Il-^ 



v^f^ = — Y 

2/^ n=0 n'=0 



n / n,n' \ rn 



(2) 



Here, the a and {3 denote sets of quantum numbers including the values of the orbital momenta 
la and and 

(3) 



" r(n + 2L + l)!^^"' 

are the functions bi-orthogonal to the Laguerre basis functions (p^f^ ) R] in the channel a 



[a) 



la+l ^-aar / 2 jJ2l a+l^ 



a^r) 



(4) 



I.e. 



= 6nn'- Here, is a scaling parameter, /i is the reduced mass of the system. 

We confine ourselves to the elastic scattering case, in which the energy e = k^ = 2fi/h^Eis the 
same in the two channels. 

The elements Uaf^ik, r) of the 2x2 matrix as the solution of the Schrodinger equation are 
represented in the J-matrix method in the form of the expansion 



(6) 



Uo.,{k,r) = Y\^<^^)bL"'\k), (5) 

n=0 

whose coefficients 6^"^^ obey the relations ||12|| : 

\ Sir' + Ci" [k) h\f, n > iV„ - 1. 

Here, are the elements of the i^'-matrix, P^n' elements of the matrix 

P = (Re - h)~^ , 

where h is the matrix of the Hamiltonian h = li^^ + li^^ + '^^V calculated in the basis 



(7) 



||q;) 0^"-'^, n = 0, Na — 1, a = 1, 2|. R is the block-diagonal overlap matrix involving the 
symmetric tridiagonal submatrices R*-"^ |15|: 



^ ^ 1 (n + 2/^ + 2)! 



(8) 



The matrix hg"'' of free Coulomb Hamiltonian 

j2 



^ la{la + 1) 2rfc 

9 I i 



(9) 



where r = Ji ■, is also symmetric and tri diagonal |T5[| : 



9la D(a) I 2rA- 



(n + 2/^ + 1)! 



n, n+1 



n+1, n 



n ^he elements of so called J— matrix 



d(") 
^ ^-n, n+1- 



2 D(a) 



(10) 



In what follows, h denotes an A/" x AT-matrix (A/" = A?"! + A^2) whose elements hij are defined 
in accordance with the rule hij = /i"^,, if z = + J2 ^a' + 1 and j = n' + J2 ^a' + 1- 

a'<a a'</3 

With the spectral decomposition of the matrix h' = ZAZ* (where superscript * denotes 
transposition) of the Hamiltonian h calculated in the orthonormalized basis ip^^'^ , n = 

0, Na — 1, a = 1, 2|, where Z is the orthogonal matrix of eigenvectors, and A is the diagonal 
matrix of the eigenvalues {Xj} of the matrix h', expression ([7p can be written as []T^ 

P = D*Z[Ie- A]"^Z*D. (12) 

Here, I is the identity matrix, D is the block-diagonal matrix involving the Na x Na submatrices 
0*-"^ a = 1, 2 of the transformation to the orthonormal basis | 'ipli'''/ , n = 0, Na ~ l}) i- e. 



(a) 



-^n+l,j 



(13) 



and 5"^°^ coincide respectively with the real and image part of the solution C^^^~^^ of 
the discrete analogue of the Coulomb Schroodinger equation M 



1 (ec.-7r/2)r -l(n+l)( 



n\ e 



r(n + /^ + 2 + ir) (2 sin^„ 



(14) 



X2-Fi(-/a + ir, n + 1; n + /c 



it: e' 



-2Wa 



Here, cTq = argT^ir + + 1), e^^" = {iqa — |)/(iga + |), = k/aa- For Z = hypergeometric 
function 2-F1 in Eq. ([l^) degenerates into a polynomial of degree /„ in e"^^^". 

The symmetric 2x2 i^-matrix corresponding to the potential (|lj), (0) has the form ||15|| : 



K 



where 



[C^_l - PJC^v].,^ = C^Nl-lik) Sa,p - Va,M JN^-l,N,ik) C'^',{k) 



Va,i3{t) = Pn^-i Na-ii^) the elements of so called P-matrix. 



mi 



(15) 



(16) 



3 Description of the method 

By choosing the orthonormal basis functions | V'i"^^! (0) in the form 

K'")-( (n + 2l + 2). ) ("-0-'e-/=L^«K.). (17) 
which is to say that the matrices D'-"-' are defined by the expressions 

Df^) = I fS-i' \ ^ (18) 



,0, i < J, 
where 

^, , (19) 

l(n + 2/„ + 2)!i ' ^ ^ 



we derive the following expressions for the P-matrix elements 



^.2(e)-(^S-.)\|^, (20) 



Pi2(e)=P2i(e) = 4^irfS-i S 



(1) ^(2) ^ ^jVi..7^Ar,.7 



It should be noted that the elements Va^pi^) (^Of ) decay (~ e^^) at the energies e > Xj\f, and 
the same is true for the i^- matrix elements (|15|). On this basis, we set ourselves the task of the 
construction of the potential (0) descriptive of the i^'- matrix in an energy interval [0, eo] 
and the system bound spectrum. In the case of arbitrary i^-matrix behaviour it is suggested 
that eo < Xj\f. 

The expressions (pOD for the P-matrix are similar in structure to those for the i?-matrix ||16 



This analogue have inspired us to implement the technique used within the finite-difference 
analogue of the i?-matrix scattering theory |13| for the retrieval of information requisite for the 
potential constructing from the ii'- matrix. For this purpose we fix some values of a^, and 
define functions 

-Puie) = — -4t^ {k) {(5«_i(fc) + C^^Umu) (5g(fc) + CP^{k)K,,) 

(21) 

V22{e) = — (A:) {{S^M.Ak) + C'^U{k)K,,) {si^l{k) + dii{k)K,,) 

A(fcj JjV2-l, Af2 

~^Nl~i{k)C\^l{k)Kl.2\ 7 

(22) 



Pi2(e) = P2i(e) = (,) , (23) 



where 

m = [S'Mlik) + C^NmKii) [S'^m + C^;v2^(A:)i^22) - C\^i{k)C'^i{k)Kl,, (24) 



by inverting expression ([15D relative to the "P-matrix elements. We set, in view of Eq. 



A, 



(25) 



where ^Xj < eo, j = 1, p| are the "internal" poles of the functions Va,f^ ( |2T| ^ 

An search for optimal combinations Cq,, A^^ is an important constituent of the method. 
Assuming that the absolute values of the off-diagonal elements of the A'-matrix are negligibly 
small, this problem can be solved in every channel a = 1, 2 separately. For this purpose, on 
every step of variational procedure involving ||^ in which the rational approximant 



^ ^ j=i e- A} 

on the interval [0, eo] of the function 



(26) 



(27) 



C^°) (fc) = cos 5, ( A;) + sin 6^ C^^^ (k) 

is constructed. Going to the next step, we increase and decrease Na. The condition that 
remainder of the approximation be small is criteria for going to new combination Oq, A^^. 

A'q, is bounded below by the requirement that the residues of Va be positive, or by the 
equivalent condition that the roots of the denominators of the functions Va alternate those of 
the numerators ||^. It is evident from asymptotic expressions for the solutions Sl°'\ as 
n 



TX L 



oo {k is bounded) |T5|: 

StKk) = sin|(T« - (n + + 1)0^ - r ln(2 n sin ^«) + ^ 

CW(^) = cosja^ -{n + la + 1)6^ - r ln(2nsin^„) + ^| 

(7^") (A;) = sin |(t„ - (n + + 1)6^ - rln(2nsin0„) + ^ + 5, 



(28) 



that if Na is large enough to fulfill the condition 
d 



d k 



Wa - {Na + L)Oa ' r\n[2 {Na - 1) sin^J} 



> 



d6n 



d k 



(29) 



the roots separation holds. 

The parameters Oq,, Na obtained above are used in the coupled-channel case. Once the 
"internal" parameters has been fixed (pSD, r = Af — p "external" triplets {Xj >eo, Z^^j, 
^JV,j, j = p + 1, A/" I remain unknown. Since the rows {Z^-^ j} and {Z^ j} of the orthogonal 
matrix Z must meet the orthonormalization conditions 

Af M 



Af 

-, '^Zn^jZj^j 



0, 



(30) 



it is evident that the minimal r is equal to 2. 

The values of the "external" parameters {Xj}, as a rule, lie in the energy range that is far 
removed from the right limit of the interval [0, eo]. If the assumption is made that for this 
energies the phase shifts and the mixing parameter are negligible, for the parameters {Xj > eo} 
the "external" parameters obtained in the construction of the rank potential in every 

channel a = 1, 2 can be used. 

There is for every bound state with the energy = — a corresponding parameter X^: 

< X^ < Xi, which fits the equation |T^]: 



det 



where 



(31) 



c(+) 



N-l 



a,l3 



(32) 



Eq. (^) is similar to the equation [TT 



det[I-G(°)(+)(iK,)V] = 0, 



(33) 



where 



a,l3\ 



G(o)(+)| 



IK,. 



9n' ) 



(34) 



are the matrix elements of the Coulomb Green's operator. In spite of the fact that Eq. (plf) 
possesses the redundant roots n = aa/'i. of multiplicity A^^^, it is preferable to Eq. (^), since 
the desired potential matrix V not enters explicitly into Eq. (^TJ). 



Phase-equivalent transformation 



From Eq. (pOD follows the expression that relates the matrices of the phase-equivalent potentials 
(for comparison see |ll4| , ^) corresponding to fixed values of a^, Na- 



V(g) = D-iQZoA(d-iQZo) -ho. (35) 

Here, A and Zo are fixed diagonal and orthogonal M x AT-matrices respectively; Q is arbitrary 
orthogonal matrix of the form 

(A'^i-th column) (A/'-th column) 



Q 




Vo ... 



(A^i-th row) 



(36) 



i; 



(A/'-th row) 



ho is the block-diagonal matrix of the free Hamiltonian involving submatrices ho""* (|T0]). Thus, 
the eigenvalues {Xj} and rows {Z^-^^}, {Zj^j} of the eigenvectors matrix Z of the truncated 



Hamiltonian matrix h' specify the potential (|l]), @ V to within the orthogonal matrix Z 
construction method from given iVi-th and A/'-th rows. Once this method has been fixed 
(e.g. 10), and the corresponding matrix has been denoted by Zq, the uncertainty of sought-for 
potential can be concentrated in the arbitrary orthogonal matrix Q (|36|). 

It should be noted that from (pO[) it also follows that the P-matrix is invariant under 
exchange of Z^^j and Zj^j signs. 



Taking into account the bound state wave function 

The components Ua of the bound state wave function are represented in the form of expansion 



Ur 



n=0 



where the coefficients for > A^^^ — 1 are given by |T2 



(37) 



(38) 



Notice that the parameters |Aj, Znij, Zj^j, j = 1, A^j determine the asymptotic properties 
of the bound state wave function. Really, from the equation ||18 



it follows that 



C«i+)(i«:) - Pn(^) J'^UN^i^^) C\^r'i^^) 



s = o, 



(1) 



^(1)(+)/ 



r] = S2/S1 



'Pui^) JnI-1,N2 



{2)(+), 
N2 ' 



To derive the expression for Si, let us define two auxiliary A/'-component vectors: 



Si d^jli 



j = 1, Ni- 1, xtvi 



.(!){+), 



d 



01 aj_^^_i 



J = A^i + 1, AA- 1, 



d' 



N2-1 



and 



9j 



With the help of (|4T|), (^2]), we can rewrite (|), (^) as 

x = Zg. 

Since Z is the orthogonal matrix, the vector x norm is equal to that of g: 



(39) 



(40) 



(41) 



(42) 



(43) 



X X 



(44) 



Notice that the functions Ua (|37|) can be rewritten in the form: 

Ml {iK, r) =Si< Xn+i ipi^> ) +Ui {in, r) ^ , 



where 



^^-2cM(+)(i^)-ca+)(i^) 



oo 
n=0 



n=0 



(a) 



oo 

E 

n=0 



4") 



Substitution of Eq. 



61) into the unit norm of the wave function condition 



oo oo 

j |ui(i/t, r)|^ dr + J \u2{ 



IK, r)\ dr 



(45) 



(46) 



(47) 



gives: 



5r 



g g 



d' 



(2)(+), 

N2-1 ' 



IK 



127 
%2-l 



+ 



/ \ui{iK,, r)\ dr + J \riU2{iK,, r)\ dr. 




(48) 



In the framework of our method the phase-equivalent transformation (^), (|36|) is a tool for 
the influence on the potential (|I]), (^ off-shell behaviour. If we only be interested in describing 
the bound state wave function properties, we may vary the first — 1 coefficients of the 
expansions ( pT] ) in view of the condition (p^. For A/" > 4, the resulting set |&^"\ n = 0, Na — 2| 

does not uniquely specify (see (ffSf)) the matrix Q (p6|), since its unknown sub matrix Q 



Q 




III 



IV 



(49) 



is of order M — 2. 



4 Results 

As an application, we invert the NN elastic scattering data. A set of Nijmegen phase shifts 
in the — 350 MeV energy range is taken as input to calculate the potentials. The obtained 
pp potentials parameters in the channels ^5*0 and F2 are listed in Table I. In Fig. 1 we 



show the phase shifts resulting from the Nijmegen partial wave analysis ^ (crosses) and the 
phase shifts obtained from our potentials (solid line). 

In the ^5*0 channel for a chosen a and N the pair of unknown "external" parameters A4, Z4 4 
has been obtained. Since the ^4^4 value is fixed by the normalization of the vector {Z^j} 
condition, only free parameter A4 remains. 

In the channel F2 the absolute values of the phase shift 63 are very small. For the 

relative deviations of the S3 values from ones of to be comparable to those for the phase shift 
61, the value of N2 must be chosen sufficiently large (see Table I). For a chosen combination of the 
tta, Ntta, we obtain three triplets of the "external" parameters {Xj, Z2J, Zgj, j = 6, 7, 8}. The 
values of Xq and A7 are obtained by making it equal to A2^'* and X^^^ calculated in the potentials 
construction in the uncoupled channels ^P2 and respectively. The parameters As, -2^2,8? -^8,8 
values were calculated in a similar way to those of the "internal"; to do this the phase shifts 
are extrapolated with functions which decay asymptotically. Finally, the orthonormalization 
conditions (pOD allow the free parameters number to be reduced to one. For the latter we used 

■2^2,7- 

The obtained values for the np potential parameters in the channel ^Si Di are given 
in Table II. The phase shifts and the mixing parameter corresponding to our potential are 
represented by solid lines in Fig. 2. 

Here, we also have obtained three triplets {Xj , Z2J, Zgj, j = 6, 7, 8} of the "external" 
parameters. The values Xq and A7 are equal to A4^'* and A4^'' obtained in the uncoupled channels 
^Di and ^5*1 respectively. Next, orthonormalization conditions (|30|) have reduced the unknown 
parameters number to four, three of which Z4 7, Z4 §, Zg § are variable. Further, the combi- 
nation of Z^^s, Zg 8 has an effect on the phase shift 60 and mixing parameter ei behaviour 
in the neighborhood of e = 0, whereas Z4 7 is responsible for ei curve form at the right of 
the interval [0, eo]. Then the parameter Ag fits Eq. ( pi]) with the deuteron binding energy 
e = -k'^ = -.0536401 /m'^. 

The deuteron wave function of the potential Nijm II is employed as input in the np 
potential constructing. It is seen from (|45D, (^) and the relation u^g\iK, r) — > e~'^^ that the 

r— >oo 

asymptotic S-state normalization As is equal to 5*1, i. e. depends on the module of the vector g 
( ^^D (see (|I8|)). The last component gg containing the parameters Z4 g, Zg g, Ag makes a major 
contribution to the vector g norm (see Table III). As is found to be the most sensitive to the 
change in Z^^g, whereas the D/S-state ratio rj [r] = —rj (|3), since u^^(i/€, r) — > — e"**^) is 

r— >oo 

strongly dependent on the value of Zg^s- 

In order to get the right behaviour of the wave function in the interior domain we use 
the phase-equivalent transformation (|35|) , (|36|). In so doing the vector x (^I]) components 
{xi, X2, X3} and {x5, Xq, x-^} (which within the factor Si are equal to the first three coefficients 
of the 5* and D waves expansions (^) on the basis orthonormal functions {\L7\j ) are varied in 
view of Eq. (P^. In Fig. 3 the deuteron wave function S and D components (functions us and 
Uo respectively) of our potential with the parameters in Tables II, III are drawn (dotted line). 
There we compare them with the wave functions of the Nijm II potential ^ (solid line). The 
deuteron parameters calculated with our wave function and those from the Nijm II potential 
are listed in Table IV. 



5 Summary 



Approximate inverse scattering method within J-matrix approach 0, ^ has been extended to 
the case of the Laguerre form factors. As a resuh the method has become apphcable to the inver- 
sion of charge particles scattering. The P-matrix has a dominant role in our procedure. Similar 
to the i?-matrix, the P-matrix contains requisite information for the potential constructing and 
on the other hand is uniquely determined by the S'-matrix. Besides, the P-matrix determines 
the asymptotic normalization constants of the bound state wave functions. 

In this work the phase-equivalent transformation fl^, ^ has been modified to be valid for 
separable potentials with the Laguerre form factors. The phase-equivalent transformation was 
used for the description of the deuteron wave function of the realistic Nijm II nucleon- 
nucleon potential . In doing this, we were interested in the behaviour of the wave function 
in the internal domain rather than in fitting the values of the quadrupole moment Qd or the 
D-state probability p^. The discrepancy between the deuteron wave functions of our potential 
and that of Nijm II potential can be removed by increasing of desired potential (0) rank 
Af. 

Acknowledgment 

It is pleasure to thank Prof. A .M. Shirokov and Prof. A. I. Mazur for helpful discussions. This 
work was supported in part by the State Program "Universities of Russia", project No 992306. 

References 

[1] N. H. Kwong, H. S. K5hler, Phys. Rev. C 55, 1650 (1997). 
[2] P. Navratil, B. R. Barrett, Phys.Rev. C 54 2986 (1996). 

[3] V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen, and J. J. de Swart, Phys. Rev. C, 
49, 2950 (1994). 

[4] Z. Papp and W. Plessas, Phys. Rev. C 54, 50 (1996). 

[5] S. A. Zaitsev, Teoret, Mat. Fiz. 115, 263 (1998) [Theor. Math. Phys. 115, 575 (1998)]. 
[6] S. A. Zaitsev, Teoret, Mat. Fiz. 121, 424 (1999) [Theor. Math. Phys. 121, 1617 (1999)]. 
[7] H. A. Yamani, L. Fishman, J. Math. Phys. 16, 410 (1975). 

[8] G. F. Filippov, I. P. Okhrimenko, Yad. Fiz. 32, 932 (1980) [Sov. J. Nucl. Phys., 32, 480 
(1980)]. 

[9] Yu. I. Nechaev and Yu. F. Smirnov, Yad. Fiz. 35, 1385 (1982) [Sov. J. Nucl. Phys., 35, 
808 (1982)]. 

[10] J. Revai, M. Sotona, J. Zofka, J.Phys. G 11, 745 (1985). 

[11] Z. Papp, J. Phys. A 20, 153 (1987); Phys. Rev. C 38, 2457 (1988). 

[12] J. T. Broad, W. P. Reinhardt, J.Phys. B 9, 1491 1976. 



[13] B. N. Zakharicv, A. A. Suzko, Direct and inverse problems. In: Potentials in quantum 
scattering. 2-nd ed. Berlin, Heidelberg, New York: Springer- Verlag, 1990. 



[14] Yu. A. Lurie and A. M. Shirokov, Izv. RAN, Ser. Fiz. 61, 2121 (1997) [Bui. Rus. Acad. 
Sci., Phys. Ser. 61, 1665 (1997)]. 

[15] J. T. Broad, Phys. Rev. A 31, 1494 (1985). 

[16] A. M. Lane, A. M. Thomas, Rev. Mod. Phys. 30, 257 (1958). 

[17] V. A. Skorobogatko, The Theory of Branching Continued Fractions and its Applications 
to Computional Mathematics [in Russian], Moscow (1983). 

[18] A. M. Shirokov, Y. F. Smirnov, S. A. Zaytsev, Teoret, Mat. Fiz. 117, 227 (1998) [Theor. 
Math. Phys. 117, 1291 (1998)]. 



Tables 

Table I. pp potentials parameters values in the channels ^5*0 and 



1^0^ a ^ 2.5fm-\ N^A 



= 1, ai = 5/m-\ Ni = 2, 
k = 3, 02 = Sfm~'^, N2 = 6 



3 






ZniJ 




A, 














1 


.1558523299 


.08517568687 


.002378372193 


-.2021238864 


.3171350908 


2 


.4062499094 


.6366888556 


.006334410556 


.2927989544 


.7864286930 


3 


.6517820375 


2.999305831 


.02352440740 


-.3709560916 


1.607103324 


4 


.6211692507 


22. 


.4982479821 


.03288283700 


2.259459783 


5 






.02790963713 


-.4402591159 


3.146276272 


6 






.01090380937 


.4980705527 


6.48 


7 






.8657366805 


.01462158020 


13.25 


8 






.02742131063 


-.5409333126 


16.20600628 



Table II. np potential parameters values in the channel ^5*1 — ^ Di. 





li = 0, Oi = 


= 2.Sfm-\ iVi 


= 4, 




Z2 = 2, 02 


= 3fm-\ N2 


= 4 


3 


ZNi,j 


Zmj 


A. 










1 


.02017320192 


2989525605 


.4678515858 


2 


.3982489552 


02212234580 


.6496845646 


3 


.01890700216 


4605017812 


1.501863086 


4 


.6535995282 


09501129016 


3.492094508 


5 


.1314698511 


5784348737 


4.238758885 


6 


.08780677309 


5929301489 


12.75 


7 


.6152722974 


05339635194 


35. 


8 


.09943159709 


006685636351 


-.03571632149 



Table III. The components of the vectors g (^2]) and x (|4l| ) in frn^^'^. 



J 


9j 


Xj 


1 


-.002345654140 


.6665 


2 


-.1089332559 


-.69 


3 


.0002896455277 


.341506600 


4 


-.03560896022 


-.1757029190 


5 


-.004683742390 


.265 


6 


-.001726119865 


.025 


7 


-.003354599822 


.045 


8 


-1.067801341 


-.01049929153 



Table IV. Deuteron properties: D/S-iatio r], asymptotic S'-state normalization As 
m~^/^, D-state probability in %, and quadrupole moment Qd in fm^. 



V As Pd Qd 



(this article) .0252 .8845 5.729 .28185 
Nijm II .0252 .8845 5.635 .2707 



Figure captions 



Fig. 1. pp phase shifts in the channels ^5*0 and — ^ F2 from the Nijmegen partial-wave 
analysis (crosses) and the potential (|I|), (0) with the parameters listed in Table I (solid line). 

Fig. 2. np phase shifts in the channel ^5*1 Di from the Nijmegen partial- wave analysis 
(crosses) and the potential @) with the parameters listed in Table II (solid line). 

Fig. 3. The S and D components of the deuteron wave functions of the Nijm II potential 
(solid line) and of the potential (P, @) with the parameters listed in Tables II, III (dotted 
line) . 
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Figure 1. 
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Figure 3. 



